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Abstract. We discuss the integrability of the Bosonic and Grassmannian massive 
Thirring models in the presence of defects through the inverse scattering approach. We 
present a general method to compute the generating functions of modified conserved 
quantities for any integrable field equation associated to the m x m spectral linear 
problem. We apply the method to derive in particular the defect contributions for the 
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^ ! 1. Introduction 

There has been a great deal of progress in studying two-dimensional integrable field 
theories in the presence of defects or internal boundaries, in both classical and quantum 
context [1, 2]. From the classical point of view, the Lagrangian formalism [3, 4] has 
shown to be significantly successful. In this framework, the usual variational principle 
^ ■ from a local Lagrangian density located at some fixed point, reveals frozen Backlund 

transformations as the defect conditions for the fields. There are many interesting 
features of these defects. It turns out that these kind of defect conditions allow for several 
types of integrable field theories [3, 4, 5], not only the energy conservation but also the 
conservation of a modified momentum, which includes a defect contribution. Moreover, 
their integrability is provided by the existence of a modified Lax pair involving a limit 
procedure, but in general it was only checked explicitly for a few conserved charges. 
As a novel feature of most of these models is that only physical fields were present 
within the formulation and therefore the associated Backlund transformations were 
called type-I [6] . However, it was noticed that not all the possible relativistic integrable 
models could be accommodate within this framework and then it was proposed a 
generalization by allowing a defect to have its own degree of freedom, and the associated 
Backlund transformations were named type-II [6]. Many examples were also discussed 
in [6] like sine/sinh-Gordon, Liouville, massive free field and Tzitzeica-Bullogh-Dodd 
model. Concerning type-II defects, it is interesting to note that for the supersymmetric 
extensions of sine-Gordon model [7, 8] and for the massive Thirring models [9, 10] those 
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auxiliary boundary fields, which correspond to the degree of freedom of the defect itself, 
had already appeared naturally. 

On the other hand, recently it was suggested an alternative and systematic new 
approach to defects in classical integrable field theories [11]. The inverse scattering 
method formalism is used and the defect conditions corresponding to frozen Backlund 
transformations are encoded in a well-known defect matrix. This matrix provided 
an elegant way to compute the modified conserved quantities, ensuring integrability. 
Using this framework the generating function for the modified conserved charges for 
any integrable evolution equation of the AKNS scheme were computed, and the type-II 
Backlund transformations for the sine-Gordon and Tzitzeica-Bullough-Dodd models has 
been also recovered [12]. It is worth noting that this method for constructing integrable 
initial boundary value problems based on the Backlund transformations has already 
been used in [13, 14]. 

The aim of this paper is to provide an alternative approach in order to establish 
the integrability of the Bosonic and Grassmannian Thirring models with type-II defects, 
which were suggested by previous approaches [9, 10]. At this point it is worth mentioning 
that there exists no totally clear definition of what integrability is for classical systems 
with infinitely many degrees of freedom. However, in this work we adopt the popular 
point of view where a system is regarded as integrable, if for its describing equations of 
motion it is possible to determine a constructive way of finding solutions and to show 
the existence of sufficient number of conserved quantities. Even though this viewpoint 
is not complete, it is sufficient for our immediate purposes. Since soliton solutions for 
these models have already been studied in [10], herein we will focus on the explicit 
construction of their conserved quantities. 

The paper is organized as follows. In section 2, for the sake of clearness we present 
the standard setting of the Lax pair approach for the general m x m spectral linear 
problem (see for example [20]). One of the most important results of this approach and 
that is the main point for oiu^ purposes is the identification of sets of coupled Riccati 
equations in order to construct conservation laws by a simple generalization of the 
method used by Wadatti and Sogo for the particular case m = 2 [15]. Other different 
generalizations, considering particularly matrix Riccati-type equations, have also been 
appeared in the literature concerning to multicomponent systems (sec for example [16] 
and references therein). 

Within this framework, for the simplest case m — 2 the defect contributions 
to the infinite set of modified integrals of motion were explicitly derived for several 
integrable equations associated with the sl{2) Lie algebra- valued Lax pair like: (m)KdV, 
NLS, Liouville and sine/ sinh- Gordon equations for which previous results derived from 
Lagrangian principles [3, 5] were also recovered. More recently, the modified energy and 
momentum for the Tzitzeica model with type-11 defects, which can be described by a 
A2 Lie algebra-valued Lax pair, were also explicitly computed [12] and showed complete 
agreement with the results obtained by Lagrangian approach [6]. Then, motivated by 
these results, in section 3 we derive a general formula to compute the defect contributions 
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to the infinite sets of modified conserved quantities for any m x m linear problem. 

In section 4, firstly we use this formula to compute explicitly the defect contributions 
to the modified energy and momentum for the Bosonic Thirring model. These significant 
results shed an interesting new light on the question of what type of Backlund 
transformations may be used as defect conditions since its Lagrangian description 
has not been derived until now. Subsequently, in section 5 we extend the procedure 
to include the pure fermionic version of the Thirring model which is described by 
5/(2, 1) algebra-valued Lax pair. We recursively compute the defect contributions to 
the modified energy and momentum and show that the results are in full agreement 
with the ones previously obtained by the Lagrangian formalism [9]. 

In section 6, we discuss in some detail how the question of involutivity of the charges 
can be addressed and how a Hamiltonian formulation can be provided to include typc-II 
Backlund transformations which carry degrees of freedom corresponding to the defect 
itself. 

Some concluding remarks are made in section 7. In Appendix A, we give a very 
short review of the sine-Gordon model with type-I defects within the framework adopted 
in this paper. 



2. Lctx pair approach 



In this section we want to discuss the main ideas of the Lax formulation in order 
to construct an infinite set of independent conserved quantities for some integrable 
evolution equations. Such equations can be formulated as a compatibility condition of 
an associated linear auxiliary problem as follows, 

(9t^(x,t; A) = V{x,t;X)^{x,t;X), (2.1a) 
aa;*(x,t;A) = ?7(x,t; A)^(x,t; A), {2.1b) 

where ^(.'z;,t;A) is an m-dimensional vector, A is a spectral parameter, and U,V are 
m X m matrices, which usually are named Lax pair. Then, from the compatibility 
condition, dxdt'^{x, t; A) = dtdx'^ix, t; A), we obtain the zero-curvature equation, 

dtU - d^V + [f/, V] = 0, (2.2) 

which gives the corresponding equations of motion for the integrable system. Now, let 
us show how to construct a generating function for the infinite set of conservation laws. 
Firstly, for every auxiliary field component with j = 1, m, we can define a set of 
(m — 1) auxiliary functions Fjj = ^j^7^ with i ^ j. Then, considering the linear system 
(2.1a,) and (2.16), it is not so difficult to identify the j-th conservation equation. 



^3 



da: 



(2.3) 



where each auxiliary functions Fjj satisfy coupled Riccati equations for the x-part, 

dx^ij — {Uij — UjjVij) -\- E^ \^ik — Tij Vkj , (2.4) 
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and respectively for the i-part, 



(2.5) 



where without loss of generality we have assumed that is a commuting field, however 
we will show later that this procedure also works in the case of anticommuting fields. 

Now, by considering solutions that vanish rapidly as |x| — )■ oo, we found that the 
corresponding j-th generating function of the conserved quantities reads, 



A wide group of integrable nonlinear evolution equations can be formulated using 
this approach, among which the most of known examples correspond to the particular 
case m = 2, e.g, the nonlinear Schrodinger equation (NLS), Kortcwcg-dc Vries (KdV) 
and the modified KdV equation (mKdV), Liouville equation, and sine/ sinh- Gordon. For 
a more complete review of these cases see for example [17]. 

It is worth noting that if the respective analytic properties of the solutions are 
considered, we can expand the functions Fjj in positive and negative powers of A and 
then solve (2.4) and (2.5) recursively for each coefficient. This immediately provides 
an expansion of the j-ih. generating function Ij in powers of A, obtaining in this way 
an infinite set of conserved quantities. In particular, the usual conserved energy and 
momentum integrals of motion, commonly also derived from the Lagrangian formalism, 
turn out to be linear combinations of these set of conserved quantities Jj, by taking into 
account coefficients for the expansions in both positive and negative powers A. However, 
these sets of conserved quantities are not functionally independent in the bulk theory 
because not all of the auxiliary fields Fj^ are independent. Although, it seems that there 
is no need to consider all the conservation laws to derive the apparently overdetermined 
sets of conserved quantities, the recent results on type-II integrable defects [12] suggest 
that for obtaining the most general form for the defect potentials, it is necessary to 
consider all the information coming from the Lax pair, i.e., all the conservation laws. 
To make it clearer, in the following section we will derive the formula for obtaining the 
modified conserved quantities which helps us to compute integrable defect potentials. 

3. Modified conserved charges from the defect matrix 

In this section we consider the Lax pair approach for constructing the infinite sets of 
modified conserved quantities in the presence of defects. 

Firstly, let us consider a defect placed at .x = 0, and suppose that there are two 
column-vector functions ^ and ^ corresponding to the auxiliary linear problems for 
X < described by the Lax pair [/, V , and for x > by ?7 and V . Then, let us 
introduce K{x,t]X), a matrix polynomial of the spectral parameter A, to connect the 
two solutions, namely. 




(2.6) 



-^{x, t; A) = K{x, t; A) *(a;, t; A) 



(3.1) 
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where, K satisfies differential equations corresponding to a gauge transformation [15] as 
follows, 

dtK ^VK- KV, d^K ^UK- UV, (3.2) 

and it is commonly named the defect matrix [3, 4, 5, 11]. This matrix is expected 
to generate the defect conditions and consequently the corresponding auto-Backlund 
transformation of the model. A classification of these defect matrices was performed as 
well as several examples corresponding to the m — 2 linear problem were examined by 
choosing a very simple form for this matrix [11]. 

We present now a straightforward extension in order to consider the mx m matrix 
auxiliary linear problem. Let us consider the generating functions (2.6) in the presence 
of the defect. 



J —ex. 



dx 



+ 



^0 



dx 



k^j 



(3.3) 



Hence, taking the time derivative and using the conservation equation (2.3), we get 

(3.4) 



d^ 
dt 



x=Q 



x=0 



Then, it is not difficult to show that the relation between the two sets of auxiliary 
functions Fj^ and Fj^ is given by 



k^j 



(3.5) 



Inserting in (3.4), one gets 



dXj Vtj 

—r- = -r-, where 
dt A. 



= Kjj + E -^■'■^ ^^J' 

k^j 



and 



3J 



33 



X V^i F,,- ] A,- + E V^iKi^ + X Vii^ik ^ky 



(3.6) 



(3.7) 



Finally, we consider the equations (2.5) and (3.2) to obtain, 

d 



dt 



Xj - In 



0, 



(3.8) 



x=0 



^33 + E ^3k^kj 

where the defect contribution to the j-th generating function of infinite conserved 
quantities is given exactly by 



^33 + E ^3k^kj 



(3.9) 



x=0 



This formula is an important result because its expansion in powers of A provides 
the defect contributions to the modified conserved quantities at all orders for every 
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mxm associated linear problem. We will apply it to study the Bosonic Thirring model, 
which is described by a sl{2) Lie algebra-valued Lax pair, and for the Grassmannian 
Thirring model, which is associated with a s/(2, 1) Lie algebra- valued Lax pair. In 
particular, it will be shown that the modified energy and momentum contributions 
can be computed from certain linear combinations of the set of conserved quantities 
Dj, taking into account all the possible conservation laws, i.e., for j = 1, m. It is 
worth noting that a similar approach was already used in [18] to prove the classical and 
quantum integrability of the sine-Gordon model with defects, by using the monodromy 
matrix language. In the same work all higher conserved quantities are found by using 
a matrix Backlund tranformation and a matrix Riccati equation. Some of these results 
are recovered using our approach in Appendix A. 

4. Case m = 2 : The Bosonic Thirring model 

In this section, we apply the method described by computing explicitly the defect 
contributions to the modified energy and momentum for the BT model. 



4.I. The bulk BT model and the linear problem 

As is well known the BT model in the bulk is integrable [19] and the associated hnear 
problem can be formulated by using 2x2 matrices valued in the sl{2) algebra as follows, 

9t*(a;, t; A) = V{x, t; A) *(a;, t; A), (4.1a) 

d^^{x, t; A) = U{x, t; A) ^{x, t; A), (4.1&) 

where the auxiliary field ^ = (\E'i, "^2)^ is a 2- vector and the Lax pair can be written 
in a compact form as, 

i [gp. - m{X' - A-2)] q{X) 

r(A) -i[gp_-rn{X'-X-')] 



U 



(4.2 a) 



(4.2&) 



-i[gp^ + m{X' + X-')] B{X) 

C(A) i [gp^ + m(A2 + A-^)] 

where for convenience we have defined p± — {(j)\(f>2^4'i(f>i), and the following functions. 



B{X) 
C{X) 



2 

iJmg 



(A01-A-V2), q{X) 
(A<^l-A-V|), r(A) 



2 

iJmg 



(A01 + A-V2), (4.3a) 
(A01 + A-VI)- (4.3&) 



2 ' ' ' 2 

Prom the zero curvature condition, we obtain the field equations for the BT model 



i{dt - d^)(j)i = m02 + g(j)l(f)2(j)i, 
i(<9i + d^)(p2 = mcpi + g(p\(j)i(j)2, 
i{dt - 5x)0l = - rn(j)l - #14^2, 

K^t + d^)(p2 = - "^01 - #2^101- 



(4.4a) 
(4.46) 
(4.4c) 
{4Ad) 
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Now, we define the auxiliary function = ^2^1 ^- Then, by using the system of hnear 
equations we have that the conservation equation can be written in the following form. 



a 



q^2i + jp- 



BT21 - —p^ 



The auxiliary function satisfies the following Riccati equations 



2 

215 



StFsi ^C+'- [gp+ + m (A^ + A'^)] T^i - BFl,. 
Now, we expand in inverse powers of A around 00, 



r2iix,t; A) = ^ 



00 -pCfc) 



21 



X,t) 



k=0 



A'^ 



(4.5) 

(4.6a) 
(4.6&) 

(4.7) 



Using the Riccati equation, each expansion coefficient r'21 {x, t) can be obtained easily 
in a recursive way. The first coefficients are given by 



-(1) 
21 



2i 



-(a,0l) +4 + ^(402)01 



m 



m 



(4.8) 



Considering, as usual, the bosonic fields (/)i{x,t) vanish at — >■ 00, the corresponding 
generating function for the conserved quantities reads 



Ii 



J — c 



dx 



q^2i + jp- 



(4.9) 



and substituting (4.7) in the expression for Ii, we get an infinite number of conserved 
quantities given by the expansion 



00 j(fc) 

k=0 

Then, the first two conserved quantities are explicitly given by 



r(0) _ ^9 



0101 + 02</'2 



? = -- r dx liMdA) - ^ (401 + 4h) 



|(0l010202) 



(4.10) 



(4.11a) 



(4.116) 



In addition, there is another set of conserved quantities that can be computed taking 
an expansion of r2i{x,t; A) in positive powers of A, 



T2,{x,t;\) = J2fi\\x,t)X' 



k=0 



In a very similar way, the first coefficients are, 

^21'* = ~ y — 025 ^21^ = 0, r^i = 4/— 

V m y m 



2i 



m 



3(a,0^)-0l-^(0l0i)0 



(4.12) 



(4.13) 
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Substituting in (4.9), we will now obtain that the conserved quantities read 

oo 



(4.14) 



A;=0 



where, the first two of them have been computed schematically, and the result is the 
following 



4 



dx 



1(2) 

^ m 



oo 

oo 



dx 



+ Y {4<i>i + 4h) + |(<^l0i402) 



(4.15 a) 



(4.156) 



Then, we have found two infinite set of independent conserved quantities as consequence 
of the two possible choices for the A-expansion of the auxiliary function r2i(a;, t; A), i.e, 
around A = and A = oo. However, these integrals of motion are not really real charges. 
Then, it is necessary to add their corresponding complex conjugate terms. In fact, these 
terms naturally raise by considering a second conservation equation that can be derived 
from the linear system (4.1a) and (4.16), namely 



r ^9 



CTi^ + —p+ 



(4.16) 



where we have introduced a new auxiliary function = ^1^2 which also satisfy a 
couple of Riccati equations. 



9.ri2 = ? + ^ [gp- - m (A' - A-')] ri2 - tV\^ , 

dtV^2 -B-'- [gp^ + m (A^ + A'^)] ri2 - CV\^. 



(4.17a) 
(4.176) 



Then, using the same scheme we can obtain recursively the first few coefficients for the 
auxiliary function ri2(a;, t;A) by considering the corresponding expansion in negative 
and positive powers of A. Doing that, the results obtained are: 

r(2) _ n r(3) _ TSi 

1) 12 — U, J- 12 — A/ — 

m \ m m 



12 



r(i) _ 
i 12 — 



r(2) 
12 



0, 



r(3) 
i 12 



m 
2i 



m 



-(5x02) - (Pl 



— (010l)02 

m 



(4.18a) 
(4.186) 



9_ 

m " ~~ \ m [m 

Prom the conservation equation (4.16), the second generating function of the conserved 
quantities can be written as follows, 

I in ' 

dx rTu-jP- ■ (4.19) 

Substituting the corresponding coefficients of the auxiliary functions for each expansion 
in A, we obtained the following conserved quantities. 



if ^ _| / dx 



15 



r(0) 
■■2 ' 



(4.20a) 
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= --/' 



dx 



-(2) 



m 



oo 

OD 



dx 



m 



401 



^l02 



) ■ 



9 

(4.20&) 
(4.20c) 



We note that the usual number of occupation, energy and momentum for the BT model 
can be expressed in the following form, 



iV = i 

25 



(if) -if) -(if) -if) 



dx 



(2) t(2)^ 



(2) t(2), 



(4.21a) 
(4.216) 



/oo r • , , 



p = 



im 
2^ 



(if)-if) + (ir-iYO 



(2) ?(2)^ 



/oo r i / 
dx - (^015x01 - 01 5a;^ 



+ 02^ 



2^x92 



4^a;02 , 



(4.21c) 



In the next section we derive the defect contribution for the energy and momentum by 
using directly the defect matrix. 



4..2. Modified integrals of motion of the BT model from the defect matrix 

As it was shown in section 3, in order to compute the defect contribution to each bulk 
integral of motion, it is necessary to know the explicit form of the elements of the defect 
matrix. Using an ansatz for the expansion of this matrix as a very simple Laurent series, 

X = i^_i + i<:o + i^i, (4.22) 

where Ki corresponds to an element of grade A', a totally consistent defect matrix was 
determined [10], which can be written in the following form 



Ae"'" -i(Aa)-^e'" 



X 



[Ae^" + i(Aa)-V"] 



(4.23) 



where 

2a = arcsin \—X^x \ , (4.24) 
-2m J 

which it turns to be related to the modified number of occupation by 4q; = gN-£,. Here, 
the boundary fields X and X'^ satisfy the following algebraic relations, 

i 



X 



He^" + 0ie- 



+ 



^1^ 



a 
1 

ia 



/)2e^" - 02e- 



4e-'" - (t>y' 



(4.25 a) 
(4.256) 



and the respective time-derivatives, 

m 
2a 



dtX 



le'" - 0ie-'« 



im 



02C'" + 02C" 
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(4.26a) 



2a 



^9 
+ 4 



(4.266) 



together with, 



m 

_ ii 

4 



im 



- 0ie 

1^1 + (t){(i)i - 4>24>2 - 4>2(f>: 



■26'" + (/)2e-'" 



) 



X, 



(4.27a) 



m 
2^ 



15' 



402 



402 



(4.276) 



where a is a real parameter. The expressions (4.25 a)-(4. 276) are the auto-Backlund 
transformations for the BT model. It was also shown in [10] that these transformations 
exhibit a complete consistency with the soliton solutions derived by applying the 
Dressing method for several transitions between them. 

Once a defect matrix is given by (4.23), the defect contribution to the modified 
conserved quantities can be calculated using (3.9). Firstly, let us consider the generating 
function of conserved quantities (4.9) in the presence of a defect, 

r-o 



/ 



dx 



q^2i + jp- 



+ 



dx 



q^2i + jP- 



Then, taking the time-derivative and using the field equations we have. 



dXi 
lit ~' 
with 



Br2i - jp+ 



BV21 - jp+ 



x=Q 



In [Kn + K,2T2i] 



x=0 



(4.28) 



(4.29) 



(4.30) 



Now, the methodology used before to obtain the conserved quantities in the bulk, can be 
applied directly. First, we calculate the conserved charges associated to the expansion 
of r2i in inverse powers of A. Then, we consider the corresponding expansion in positive 
powers of A. Implementing this procedure, we find out the following results. 



(2) 



,2ia 



a 



+ ^X(f)le 
m 



(2) 



lae 



-2ia 



lag 



m 



X0^e 



Performing the same procedure for the generating function (4.19), one gets 
dXs r~~ , if ~ ] r 1 _ dD2 

— = cri2 + — p+ - c;ri2 + — p+ = — —, 

dt 4 x=o 4 x=o dt 



{AM a) 
(4.316) 

(4.32) 
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where 



L>2= - In [K22 + i^2iri2] 



x=0 



and from which, we obtain the following coefficients, 



D. 



(2) 



e 

a 



-2ia 



+ ^XVie- 
m 



D 



(2) 



m 



'2 6" 



(4.33) 

(4.34a) 

(4.346) 



By analogy we find that defect energy and momentum for the BT model can be written 
in the following way. 



im 

Jz/r) — 

2^ 

+ 
im 

- — 

^9 



J 2g \ aj ^ 
e- - Xt^i e--) - I (x4 e-ia ^ ^t^^ e^) , 



- Df^) + (Df ^ - Df^)\ = ^ ( a - ^ ) (e^'" + e'^^") 



+ {X(t>\ e'" - Xt0i e-^") + I [x<pI e-- + Xt^^ e^") 



(4.35 a) 



(4.356) 



We can note that it is possible to rewrite these results in an alternative form by using 
the Backlund transformations (4.25a)-(4.276), as follows 

• r 

a. 



E 



D 



D 



2ia 



2ia 



+ 



m 

9 
m 



1 



a 1 cos (2a 

a 



cos {2a) , (4.36a) 
(4.366) 



These expressions for the defect energy and momentum seem not to have been reported 
elsewhere in the literature and constitute a very important result in order to address in 
future works the question of the Lagrangian formalism as well as quantum aspects like 
the transmission matrix. 

Since the intcgrablc defect conditions for the BT model have already 
been determined by giving the corresponding auto-Backlund transformations, the 
integrability of the model in the presence of defects, following the intcgrability criteria 
adopted in this work, is provided by the existence of the defect matrix and the explicit 
computations of the modified conserved quantities. 



5. Case m — 3 : The Grassmanian Thirring model 

In this section, we will apply the formalism on the GT model. As it was done in the 
Bosonic case, we will construct an infinite set of independent conserved quantities for 
the bulk theory, and then derive the corresponding defect contributions to the modified 
conserved quantities. 
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5.1. The hulk GT model and the linear problem 



The equations of motion for the GT model can be described as a compatibihty condition 
of the following associated hnear problem, 



(5.1a) 
ib.lh) 

We remark that the 



d^^!{x, t; A) = U{x, t; A)^(x, t; A), 
dt'^{x,t;X) =V{x,t;X)'^{x,t;X). 

where the three-wave function \1' has the form (^1,^2,^3) 
Bosonic or Grassmaniann property of these functions depends on the conservation law 
that will be considered. The Lax pair U, V arc 3x3 matrices belong to the s/(2, 1) Lie 
algebra (see for example Appendix G in [12]), and can be explicitly written as follows, 

f p_ + f (A2 - A-2) gi(A) 

_i^_ + i2i(A2-A-2) q,{X) , (5.2a) 

ri(A) r2(A) im{X^ - A'^) 



U 



-A+f(A2 + A-2) Si (A) 

A+f{X^ + X-^) B2{X) 

Ci(A) C2(A) im(A2 + A-2) 

where just for simplicity, we have also defined the functions 



(5.2&) 



qi= -C2 = -iy^(AV^i + A-V2), 



-Bo 



12 



r2 



(5.3) 

= i^ (a^I - A- VI) , (5.4) 
A-V2) . (5.5) 



mg 



(A^i- 



By applying the zero-curvature condition we obtain the field equations for GT model. 



i{dt - 5^)V^i = mV'2 + giplip2ipi , 
i{dt + dx)i>2 = rni^i + #IV^iV'2 , 
\{dt - 5^)V^I = - - #2V'2V^i , 

In components, the set of differential equations (5.1a) and (5.16) read. 



a.*2 = 



10 /,o , ox 



*i +91*3, 



im 



*2 + g2*3, 



5^*3 = ri*i + r2^'2 + (A^ - A ^) *3, 



and 



5**1 = - 



^ - Y (A^ + A-) 



^1 + Si^3, 



(5.6) 
(5.7) 
(5.8) 
(5.9) 

(5.10a) 

(5.105) 
(5.10c) 

(5.11a) 
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^^2 + ^2*3, 



(5.115) 



= Ci^i + C2^2 + im (A' + A"') ^s- (5.11c) 

Now, by defining the auxiliary functions r2i = ^2^r^ ^31 = ^3^]^^, we obtain a 
first conservation equation from (5.10a) and (5.11a), namely. 



a 



qi^3i + yp- 



= [B,r,, - A] , (5.12) 

where r2i and T^i satisfy the following coupled Riccati equations for the x-part, 

9x^21= - (ifl'P-)r2i + g2r3i - gir2ir3i, (5.13a) 

dxTsi = ri + rsFsi - ^ [gp- - m (A' - A'^)] r3i, (5.136) 
and for the t-part, 

dtV^i = 2Ar2i + S2r3i - Sir2ir3i, (5.i4a) 



dt^3l — Ci + C2r2i + 



. irn ,.j . 2\ 
A+ — (A^ + A ^) 



31 



(5.146) 



Now, by firstly considering an expansion in inverse powers of A for the auxiliary functions 
as. 



00 



Tij{x,t]\) = J2 



X, t) 



k=l 



A^ 



(5.15) 



and inserting this in the Riccati equations (5.13a) and (5.136) we find that the first 
coefficients of the expansion are given by 



31 



r(3) _ 

^ 31 — 




m 



2i 



.(2) 
31 



m 



1 J- 21; 



m 




where r2i'' and T21' satisfy the following differential equations, 

d^rS^ -ig[4i^2 + 4i^i)rS, 

d^rf^ ^ -ig (^|V^2 + V'lVi) rg^ + ^ [4d.4) + Si^V'iV: 



m 

(2) 



(5.16a) 
(5.166) 

(5.17a) 
(5.176) 



Then, we find out in this case that the associated generating function of the conserved 
quantities are given by. 



f 

J — ( 



dx 



9ir3i + —p- 



(5.18) 



and substituting the respective coefficients for the auxiliary function r3i in the expansion 
in A, we found that the lowest conserved quantities are given, 

r(0) _ il 
2 „ 



dx 



dx 



^1^2 + -01^1 



m 



m 



(5.19 a) 
(5.196) 
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Now, to compute a second infinite set of conserved quantities we have to expand the 
auxiliary functions around A = 0, i.e, in positive powers of A, as follows 



(5.20) 



k=l 



By following the same procedure, we obtain that the respective first few coefficients for 
each expansion are given by. 



f(i) 

^ 31 



f(3) 
31 



with 




4, 



^(2) 
31 




1p2 



r(i) 

^ 21 ) 



2i 
m 




a.fg) = ig (^1^2 + V'lVi) rg^ - ^ {4d.4) - ^^944- 



(5.21a) 
(5.216) 

(5.22a) 

(5.226) 



Then, we have that the corresponding first charges associated to this expansion of Fsi 
are given as follows. 



(0) 



r(2) 



2 

oo 



dx 



dx 



42-1^2 + 4lfpl 



(5.23a) 



(5.236) 



Clearly, these charges are not totally real and therefore it is necessary to add the 
hermitian conjugate terms. To do that, we need to consider other contributions coming 
from two more conservation equations that can be derived using (5.106), (5.10c), (5.116) 
and (5.11c), namely 



ig 

q2^32 - ^P- 



23 



da: 



2J- 32 



A 



23 



(5.24a) 
(5.246) 



where we have introduced some other auxiliary functions = ^1^2^^, = ^3^^^, 
Tis = ^'i^'s ^, and r23 = \1'2^1'3 ^. It is very easy to check that the set of Riccati equations 
satisfied by these auxiliary functions can be written as, 

dx^l2 = ifi'P-ri2 + gir32 - g2ri2r32, (5.25a) 

d^r32 = r2 + riri2 + {gp- + m(A' - A"')) Tas, 



^ris = 5i + ^ [gp- - m (A^ - A ^)] ri3 + r2ri3r23, 
5.r23 = 92 - ^ [gp- + m (A^ - A-')] r23 - riri3r23, 



(5.256) 
(5.25c) 
(5.25cO 
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dt^i2— — 2Ari2 + -Bir32 + B2ri2r32, 



dt^32 — C2 + Ciri2 



+ C2ri3r23 



ar 



't^ 23 



Bo + 



r23 — C'iri3r 



23- 
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(5.26a) 
(5.266) 

(5.26c) 

(5.26cO 



Now, these equations are solved by expanding each of the auxihary functions in positive 
and negative powers of the spectral parameter A. Performing similar computations, the 
first few coefficients for these auxiliary functions can be determined, and the results 
read 



^ 23 



^(1) 
23 





-(1) 
23 



r(2) 

^ 32 



-(1) 
32 




^1 



^(1) 
13 



_r(i) 

^ 32 





^(2) 
32 




p(2) _ p(2) _ U2) _ U2) _ „ 
^ 13 ~ ^ 23 — ^ 13 — ^ 23 — 

and 



.(3) 
32 



^(3) 
32 



.(3) 
13 



^(3) 
13 



r(3) 

^ 23 



^(3) 
23 



m 

2 

m 

2 
m 

2 
m 

2 
m 

2 
m 




—I 



{dM + (^2 + ri22) + (V1V^2) V'l 

(5.V'2) + (V'l - V'l ri22) + f (^iVl) ^2 

m ^^^^^^ 
2^ 




m 
2^ 




-u/ - (a.v^2) + + (^i^O ^2 

m v2 2\m\ / 



m V2 2\ m \ / 




2 2 V m 



(^iVi) V'l 



together with the following relations, 

a,rg = i^ (^iVi + V'IV2) rg, 
5,rS'2^= - i^ (v^.Vi + ^^^2) f S^2\ 
a.rSJ = i^ (v'lVi + t/;|V2) rSJ + ^ (V'la.Vi) - 2i^ ii^M 



a.fSJ = - i^ (v'lVi + 4^2) r 



m 

(2) _ 2^ 
^2 m 



{lp2dxtp2) + 215- (lAl^2) 



(5.27a) 

(5.276) 
(5.27c) 

(5.28a) 
(5.286) 
(5.28c) 
(5.28(/) 
(5.28e) 
(5.28i) 

(5.29a) 
(5.296) 
(5.29c) 
(5.29(i) 
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Now, wc will compute the corresponding conserved quantities from the conservation 
equations (5.24a) and (5.246), namely 



dx 



oo 
oo 



32 



-P- 



dx 



23 



(5.30a) 



(5.306) 



Therefore, by a straightforward substitution of the each expansion coefficient, we easily 
get the following results. 



r(0) 



2 i-< 



r(2) 



dx 



r(2) 



-L 



oo 
oo 



r(2) 



r(2) 



dx 



dx 



-oo 
oo 



dx 



dx 

m 

_^ 
m 

m 
2<7 



- I 



(0) 
2 ' 



r(0) 







f(o) 

■■3 ' 



+ 



15 ( V'lV'l + '0lV'2 



m 



V^^V'2V'iV'i) 

V'2^2V'l'i/'l) 



(V-I^xV-i + V'l^xV'l) + 2i5 (^IV'i + ^-1^2) + ^ (4V'2V'IV'i) 



m 
2i5^ 



(5.31a) 
(5.315) 
(5.31c) 
(5.31d) 
(5.31e) 



Then, from all these conserved quantities together with the ones derived in (5.19a), 
(5.196), (5.23a), and (5.236), we can notice that 

(5.32) 



r(") 



(n) _ _(n) 



+ ir = I 



r(") 



3 ' 



in) 



Therefore, it is convenient to define the following quantities, 
l(o) = (lf)-lW_i(o))^ J(0) ._ 



j(0) jiO) ^(0)^ 
-■■1 ^2 -"-a y 



I(^) = (li^) + if + if), 1(2) = (l(^) + if + IP)^ 



(5.33a) 

1 -rx2 -rx3 ii — ^^x^ -r ±2 -^3 ;) (5.336) 

in order to get the usual conserved number of occupation, energy and momentum for 
the GT model by performing a simple combination, namely 



2ig 



E 



m 
8i^ 



1(0) _f(o) 

1(2) _f(2) 
1(2) +f(2)' 



dx 
dx 



dx 



(5.34a) 



(5.346) 
(5.34c) 



8ig 

In the next section, we show how this framework can be used in order to compute the 
modified conserved quantities from the defect matrix for the GT model. 



5.2. Modified integrals of motion from the defect matrix 

Let us implement a defect placed at the origin x = 0, and the relation between the 
respective auxiliary wavefunctions in x < 0, and the corresponding for a; > 0, is 
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given by 

t] A) = Kq{x, t; A) t; A), (5.35) 

where the defect matrix Kq satisfied the gauge transformations, 

dtKc =VKg-KgV, (5.36 a) 

d,KG = UKG-KGU. (5.365) 

The exphcit form of the defect matrix Kq was also computed in [10] and can be written 
in the following simple form, 





Kc 





m 

where the elements K± are given by, 
K± — Aexp 



2m 



i(Aa) ^ exp 



m 
Y m 

-(A + i(Aa)-i) 



2m 



(5.37) 



[A - i(Aa)-i] ± ^ [A + i(Aa)-^] X^X, 



(5.38) 



and the Grassmannian boundary fields X and X^ satisfy the following defect conditions. 



za-^(^2-^2)-7^XtX^2, 
2m 



X ={i;, + ^,) + '^ij,X^X 

2m 

together with their respective time derivatives. 



dtX = 2^(V'l - V'l) - + V^2) - J [i^Ul + 4li^l + + 4^2 



a,xt = ^{4- 4) + Y^42 + 4) + f [4ii^i + ^iVi + 42i^2 + ^IV's 

and their x-derivatives. 



dxX = — (V'l - A) + + V'2) - J [4lA + 4lA - 42'^2 - 42'^2 



X, 



d^x^ = ^(^I - 4,) -'-^{42 + 42) + y 4,^1 + ^iVi - V^lV'2 - 42i^2 x\ 

za z 4 L J 



1^ 



(5.39) 

(5.40 a) 
(5.405) 

(5.41a) 
(5.415) 



which correspond precisely to the auto-Backlund transformations for the classical GT 
model [9, 21]. 

Let us consider now the defect contributions to the conserved quantities. As we 
have discussed across this work, the entries of the defect matrix determine the modified 
conserved quantities from (3.9). First of all, let us consider the first set of conserved 
quantities given by (5.18) in the presence of a defect, 

r-o 



- 

J — < 



+ 



dx 



(5.42) 
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Then, taking the time-derivative and using the formula we found that Ii + Di is 
conserved, where the defect contribution Di to this first set of conserved quantities 
is exphcitly given by 



Di = - In 



Ku + Xi2r2i + KisFsi 



x=0 



(5.43) 



Hence, by taking the both expansions in negative and positive powers of A and the 
exphcit form of the defect matrix (5.37), we get 



(0) 



2m J 



X^X. 



(5.44a) 



= - ^X^X - ^x4, d[ 



9 



2g 



)(2) 



-Xi^l 



(5.446) 

mm mm 
In a similar way, repeating the computations for the other two generating functions 
(5.30a) and (5.306), we find that the respective defect contributions are given by. 



L>2 = - In [K2iru + K22 + i^23r32] 

D3 = - In [KsiTis + i^32r23 + i^aa] 
Using them, we obtain 



x=0 



x=Q 



D. 



(0) 



'-^) X^X, = ( '-^) X^X, = = 5 



2m J 



2m J 



(0) 

'a ) 



D. 



(2) 



m 



(2) 



m 



D?^ = - ^-ixH^ - ^x^L 



D. 



5(2) _ 



m m 
m m 



As it was expected, we also have the relations D^'^ = Dj^"-* H-Dg"'' ^-nd Dg"^ 
Then, defining by analogy the following defect quantities. 



p(0) ^ ^(0) _ 



)(0) 



= 



iga 



m 



X^X, 



p(o) _ 5(0) _ Q 



D. 



(0) 



(2) 



D 



(2) 



D 



(2) 



^ ) X^X, 
m ) 

m m 



p(2) ^ 5(2) + 5(2) + 5(2) 



4ia5f 



m m 



(5.45a) 
(5.456) 

(5.46a) 

(5.466) 
(5.46c) 

5(-)+5w. 

(5.47a) 

(5.476) 

(5.47c) 
(5.47d) 



the corresponding defect number of occupation, energy and momentum can be written 
in the following way, 

2\g \ / m 



Ed 



PO-^g (^^'^ + ^^'0 ^ ^ [^^^^^ + ^^^^^ + ^^^^^^^ " ^^^^ 



(5.48a) 
(5.486) 
(5.48 c) 
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Notice that, we can rewrite these results by using the Backlund transformations (??) 
and (5.39) in a more convenient form, by ehminating the auxihary fields X and . 
The results are the following. 



E 



D 



ag 



9 



9 



(4^24^2), (5.49) 
(V^|V^2V'|V'2). (5.50) 



Then, we have particularly derived in an alternative way the defect energy and 
momentum for the Grassmannian Thirring model in the presence of type-II defects. 
These results are in complete agreement with the ones obtained based on variational 
principles [9]. The advantage of the approach used in this work is that it provides us 
a formal way to compute explicitly an infinite number of conserved quantities ensuring 
integrability of these models. However, taking into account recent developments in the 
area of integrable defects [22] it turns out interesting to investigate how to provide a 
Hamiltonian formulation to include the models discussed in this work since the role of 
the degree of freedom corresponding to the defect itself needs to be clarified. In the next 
section, we discuss some useful ideas in addressing this question but a more complete 
description of the Hamiltonian approach will be postponed for a future work. 



6. Comments on Liouville integrability 

So far, an infinite set of independent modified conserved quantities arising from the 
defect contributions have been systematically constructed through a general formula 
derived from a variant of the classical inverse scattering method, which are from 
our point of view sufficient for these kind of defects to be regarded as integrable. 
However, the question of the involutivity of such quantities (required to discuss complete 
integrabihty in the sense of Liouville) still has to be answered. In this section we address 
some relevant facts related with this problem. 

Certainly, the Hamiltonian formulation of the classical inverse scattering method, 
which is essentially based on the concept of a classical r-matrix [23], is perhaps the most 
elegant and convenient framework to discuss involutivity. Let us start with the main 
aspects of the method in order to discuss this issue in the bulk. In the inverse scattering 
method the construction of the action-angle variables depends basically on the entries 
of the monodromy matrix t(A) — T(oo, — oo; A), where 

T{x,y;X)^rexp\ [ U{z;X)dz\, (6.1) 



y 



is the transition matrix, U{x;X) is the x-part of the Lax (2.16) at a given time, and 
V being the path ordering. As it was noticed in [23, 24] , the existence of the classical 
r-matrix, an x matrix which satisfies the relation 

{U{x- Ai) ® U{y- A2)} = 6{x - y) [r(Ai, A2), U{x; X,) ® I2 + I2 ® U{y- X,)] , (6.2) 
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permits us to write down the Poisson brackets between matrix elements of the transition 
matrix in the following form, 

{T{x, y- Ai) <f) T{x, y- M)} = [r(Ai, A2), T{x, y; Ai) ® T{x, y; A2)] , (6.3) 

from which it is derived that logarithm of the traces of the monodromy matrix commute 
for different values of the spectral parameter, namely 

{lnr(Ai),lnr(A2)} = 0. (6.4) 

Expanding (6.4) with respect to Ai and A2, wc get the involutivity of the conserved 
quantities {In}, which means that r(A) is the generating functional for the integrals of 
motion. Let us mention that the explicit form of the r-matrix for the massive Thirring 
models can be found in [23] . 

Now, let us discuss how the classical r-matrix approach is modified by including 
jump-defect (or point like-defect) in the system. As it was noticed in [18] and more 
recently in [22] , the description of an integrable defect in the r-matrix approach requires 
to introduce a modified transition matrix, 

r{x, y; A) = T{x, 0+; A) ^"^(0; A) f (0", y; A), (6.5) 

which is a combined bulk-defect transition matrix, where r(a;,0"'";A) and T{0^,y;X) 
are the bulk transition matrices corresponding to a; > and a; < respectively, and 
K{X) = K{0; A) is the defect matrix whose entries are evaluated in the single point 
X = 0. The key point in order to show Liouville integrability is to require that the 
defect matrix satisfy the Poisson algebra (6.7), namely. 



{K-\X,) (f) K-\X2)} = [r(Ai, A2), K-\Xi) (2) K-\X2) 



(6.6) 



where r(Ai,A2) is the same classical r-matrix for the bulk transition matrices. Hence, 
the above requirement is a sufficient condition to obtain the important result, 

{T {x, y; Ai) ® r {x, y; A2)} = [r(Ai, A2) , r{x, y; X,) T {x, y; A2)] , (6.7) 

which guarantees the existence of the infinite set of modified conserved quantities. 
Similar to the bulk theory, the explicit form of these integrals of motion can be extracted 
by introducing the following representation for the bulk transition matrix [25], 

T{x, y; A) = (1 + W{x; A)) e^^^'^^^) (1 + W{y; A))"^ , (6.8) 

where W{x; A) is an off-diagonal and Z{x, y; A) a diagonal matrix. Then, the logarithm 
of the trace of the modified monodromy matrix (6.5) is the generating function of 
the modified conserved quantities, where the defect contributions in an appropriate 
expansion in A, read [18, 22]: 



D{X) = In 



(1 + T^(0+, A)) ' K-\X) (1 + W{0-, A) 



(6.9) 



where the subscript ii denotes the leading term coming from the trace of the modified 
monodromy matrix for the given expansion. At first sight, it seems not to exist a direct 
relationship between the above result and the generating function (3.9) what we have 
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worked with. However, note that W{x] A) satisfy a matrix Riccati equation similar 
to (2.4), which permits us to derive recursively its coefficients in an asymptotic series 
expansion as A — )■ oo and A — )■ 0, and to demonstrate order by order that the results 
are complete equivalent for the massive Thirring models. This analysis deserves more 
attention than we could give at this moment and is left for future investigations. 

It is remarkable that the approach we have adopted (see also [12]), uses essentially 
an on-shell defect matrix which implies that its entries have non- vanishing Poisson 
brackets with the bulk monodromy matrices elements. This fact has already been 
outlined in [6] for the Hamiltonian formulation of the type-II defects in the sine- 
Gordon and Tzitzeica model, where the defect conditions appear as a set of second class 
constraints on the fields, which induces a slight modification of the canonical Poisson 
brackets. This issue indeed can be solved by working firstly with the off-shell defect 
matrix to compute the Poisson brackets and then derive the constraints as consistency 
conditions in constructing the time-like operator in the Lax pair such that the zero 
curvature condition provides the same equations of motion as the ones coming from the 
Hamiltonian evolution derived via Poisson brackets [22] . 

7. Concluding remarks and perspectives 

In this paper, we have presented a systematic approach to the integrability problem of 
the Bosonic and Grassmannian Thirring models with typc-II defects using the inverse 
scattering formalism. Such a formulation allows to compute explicitly the modified 
conserved quantities to all orders in terms of the defect matrix. 

We have followed the approach to defects in classical integrable field theories given 
in [11], to present for the mxm linear problem, a direct generalization for the generating 
function of the defect contributions for the modified conserved quantities to all orders. 
We have successfully applied this procedure in the case of Bosonic (m = 2), and 
Grassmannian (m = 3) Thirring models. For the latter, we have recovered previous 
results obtained from the Lagrangian approach [9]. For the Bosonic Thirring model, 
we have also derived explicitly the defect contributions for the energy and momentum. 
These results seem not to be reported elsewhere in the literature, and should be of 
crucial importance for further studies on its possible Lagrangian description. 

However, a remarkable aspect is that despite the duality relation between the sine- 
Gordon and Thirring models, it is not clear for the author why they allow different 
types of integrable defects, in the sense that the sine-Gordon model allows type-I and 
type-II defects, but the Thirring model apparently only allows type-II integrable defects. 
It is worth exploring in more detail the relation between these type-II defects through 
bosonization techniques in future developments. 

Finally, taking into account the reasons already mentioned in section 6 about the 
classical r-matrix description of Liouville-integrable point-like defects in integrable field 
theories, it should be interesting to investigate the aspects of the complete integrability 
of the Bosonic and Grassmannian Thirring models with type-II defects within the 
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Hamiltonian framework. Perhaps the most important motivation to perform that further 
study is the possibihty of considering the quantization by the so-called quantum inverse 
scattering method [24, 26, 27]. Some of these questions are expected to be developed in 
future investigations. 
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Appendix A. A brief review of Type-I defect sine-Gordon model 



The sine-Gordon equation 

- dl(p = -m^ sin (p, (A.l) 

can be derived as a compatibility condition for the associated linear problem given by 
the Lax pair 

r(A) iidtip) 
where it has been defined the fields. 



V 



-{{d^cp) A{X) 
BiX) iid^p) 



(A.2) 



q{X) = -_(Ae^-A-V^), 
A{X)= -^(Aef +A-Vf), 



r(A) = -(Ae-^-A-^e^), 



B{X) = ^(Ae-f + A-^e^). 



From the linear system, we can derive two conservation equations, namely 



d. 



rri2 + \ {dt^) 



- - {d,<p) 



5ri2 + - {d,ip) 



(A.3) 
(A.4) 

(A.5) 
(A.6) 



where the auxiliary functions = $2^1 ^ and ri2 = $i$2 ^ been introduced, which 
satisfy the Riccati equations. 



d,r,2 = ? - ^ i9t^) ri2 - rvl^. 



(A.7) 
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Solving these equations for and ri2, by considering the respective expansions in 
both positive and negative powers of A, we get 



^(0) 
21 

(1) 
21 

(1) 
12 

(2) 
21 

(2) 
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— [dt^ + d^(p] e 
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m 



— [dt(p + d^(p] e^, 
m 



(A.8) 
(A.9) 

(A.IO) 

(A.ll) 

(A.12) 

(A.13) 

v^t-r ^x-r I I y • (A. 14) 

Therefore, from the two generating functions of the infinite conserved quantities 
determined by (A. 5) and (A. 6), we have that the first non- vanishing conserved quantities 
are given by 
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(A.15) 
(A.16) 
(A.17) 
(A.18) 



Therefore, the usual energy and momentum for the sine- Gordon model are recovered by 
combining these results as follows. 



= im (l5') - I^'^ + if = p + } -m'cos(/p 



dx {dt(p){d^>p). 



(A.19) 
(A.20) 



Now, considering the generating function of infinite charges in the presence of a 
defect, we have that the first contributions to the modified conserved quantities are 
given by 



x=0 



Di= - In Kii + 7^12 V21 
where the defect matrix can be explicitly written as [4], 



(A.21) 



(A.22) 
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Hence, taking into account the expansion in both negative and positive powers of A and 
the form of the defect matrix (A. 22), we found that 



(1) 



icre 



a m 



Now, following the same procedure for the second generating function, one gets 



Do^ - In 



i^2iri2 + K; 



22 



x=0 



from which we obtain the following contributions, 

_igi(^-^)/2 _ j_ 



a m 

Then, the corresponding defect energy and momentum are given by. 



En — im 



Pn — im 



(1) n(i) 



1) 



D\^' - D),' + D\'' - D),' = 2m 



(TCOS 



(TCOS 



2 

(^ + 



1 

— COS 

a V 2 



1 (^-^ 

— cos 

a V 2 



(A.23) 

(A.24) 

(A.25) 

, (A.26) 
, (A.27) 



which are in complete agreement with the results previously obtained from both the 
Lagrangian [4] and the r-matrix [18] approach. 
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